Structure of the exotic spin-flop states in BaCu 2 Si207. 
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The unusual 2-stage spin flop transition in BaCu2Si207 is studied by single-crystal neutron diffrac- 
tion. The magnetic structures of the various spin-flop phases are determined. The results appear 
to be inconsistent with the previously proposed theoretical explanation of the 2-stage transition. 



I. INTRODUCTION 



The quasi-one-dimensional (quasi-lD) S = 1/2 an- 
tiferromagnet BaCu 2 Si207 is recognized as an almost 
ideal model material for studying exotic spin dynamics 
in weakJUfpinteracting quantum spin chains.EJ In a series of 
papersatffl we have demonstrated that in this compound 
it is possible to observe a separation between single- 
particle spin wave excitations, characteristic of a classical 
system with long-range magnetic order, and a gapped 
continuum of states that is a property of the quantum 
ID Heisenberg model. Quite recently however, another 
aspect of the physics of BaCu2Si2C>7 attracted a great 
deal of attention. It was found that in the ordered state 
of BaCu 2 Si207, an application of a magnetic field along 
the direction of ordered moment leads to two consecutive 
spin re-orientation transitions a This phenomenon is in 
contrast with the typical behavior of a conventional an- 
tiferromagnet in a field, where only a single "spin-flop" 
transition occurs. 

By far, the largest energy scale of magnetic interac- 
tions in BaCu2Si2C>7 is that of antiferromagnetic nearest- 
neighbor coupling within the spin chains, that run along 
the c axis of the orthorhombic crystal structure (space 
group Prima, o- = 6.862(2) A, b = 13.178(1) A and 
c = 6.897(1) A. O). The in-chain exchange constant is 
,/ = 24.1 meV, as previously determined from inelas- 
tic neutron scattering-, experiments and bulk suscepti- 
bility measurements. lira Exchange coupling between the 
chains is much weaker, with J x = —0.46 meV and 
J y = 0.20 meV, for nearest-neighbors along the a and 
b axes, respectively, and J3 = 0.15 meV along the 
(110) diagonal.El The spin wave spectrum shows a small 
anisotropy gap of about 0.3 meV, corresponding to a: 
easy axis along the (001) crystallographic direction 
Long-range magnetic ordering occurs at a low Neel tem- 
perature of T N ss 9.2 K.u The ordered magnetic moment 
is rather small, mo =0.15 (13 per Cu 2+ .BEl In the or- 
dered state the spins are parallel to the magnetic easy 
axis. Nearest-neighbor spins along the b and c axes are 
aligned antiparallel relative to each other, while .nearest 
neighbors along a are ferromagnetically aligned Bu The 




FIG. 1: A schematic view of the zero-field magnetic structure 
of BaCu2Si207. The S — 1/2 Cu 2+ ions form slightly zigzag 
antiferromagnetic chains. 



spin structure is visualized in Fig. [l]. The phase transi- 
tions that are the subject of this work occur at relatively 
low fields H cl = 2.0 T and H c2 = 4.9 T, applied along 
the easy axis, and produce characteristic jumps in the 
magnetization curve.El To date, the spin structure in the 
high-field phases has not been determined. However, in 
Ref. H it was suggested that upon going through the first 
transition, the spins "flop" into the (a, b) crystallographic 
plane. The second transition was then described as an 
unusual spin rotation within the plane perpendicular to 
the field, and is somewhat similar to that recently found 
in K 2 V 3 8 £ 

Based on an analysis, of the high-field magnetization 
data, it was conjecturedQ that the unusual second spin re- 
orientation in BaCu2Si2C>7 is a result of competition be- 
tween Dzyaloshinskii-Moriya (DM) off-diagonal exchange 
interaction in the spin chains and isotropic exchange cou- 
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pling between the chains. In the present work we re- 
port the results of magnetic neutron diffraction studies 
of the various spin-flop states in BaCu2Si207. We find 
the results to be inconsistent with the previously proposed 
model. It appears that inter-chain coupling dominates 
over DM interactions, and that the series of spin re- 
orientation transitions can not result from a competition 
between the two effects alone. Alternative mechanisms 
for the two-stage spin flop transition in BaCu2Si207 are 
discussed. 



II. EXPERIMENTAL 

A 5 x 5 x 4 mm 3 sample was mounted in a 6 T cryomag- 
net on the D23 lifting counter diffractometer at Institut 
Laue Langevin. The (0,0,1) axis was aligned parallel 
to the field direction. With 8 magnetic Cu 2+ per unit 
cell, in the ordered state all magnetic reflections coincide 
with nuclear ones. In addition, the ordered magnetic 
moment in BaCu2Si20y is small, and detecting the mag- 
netic contributions to Bragg intensities is rather difficult. 
The measurements were therefore performed in differen- 
tial mode. At H = T, 4 T (below H c2 ) and 5 T (above 
H C 2) the diffractometer was consecutively positioned on 
each accessible Bragg reflection. Without moving any in- 
strument motors, the peak intensity was measured first 
at 2 K and then again at 10 K. This measurement yielded 
the ratio of nuclear and magnetic scattering intensities. 
The absolute values of nuclear intensities were separately 
measured at 10 K in standard rocking scans. The rock- 
ing curves also provided an estimate of the background 
around each Bragg position. The measured integrated 
nuclear intensities, backgrounds, and peak intensity ra- 
tios were then used to estimate integrated magnetic in- 
tensities. In another experiment, the peak intensities of 
(411), (031) and (051) reflections were monitored as a 
function of applied field in the range 0-5.5 T. Peak in- 
tensities measured in this fashion were assumed to be 
proportional to integrated intensities, with the normal- 
ization factor determined at H = 0, as described above. 



Phase II 
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Phase III 



FIG. 2: A model for the spin arrangement in BaCu2Si207 in 
magnetic fields applied along the c axis: Phase II for 2 T = 
H cl < H < H c2 = 4.7 T and Phase III for H > H c2 . The 
spins are largely confined to the (a, b) crystallographic plane. 



at each field was reduced to determining the orientations 
of four spins, one from each chain, labeled as S1-S4 in 
Fig. @. 



A. Spin structure at zero applied field: Phase I 

As a first step, we have verified that our new data 
collected at H = are consistent with the structure 
previously reported in Ref. |[ The spins are parallel to 
the c-axis, with nearest neighbor spins along the c and 
b axes antiparallel to each other, and a parallel align- 
ment of nearest neighbor spins along the a axis, as shown 
111 Fig. 0. This type of spin arrangement minimizes all 
inter-chain interaction energies that were independently 
measured using inelastic neutron scattering] 3 ] Magnetic 
Bragg intensities calculated for this c-axis colinear state 
are in good agreement with those measured at H = in 
the present work. This comparison is made in Table 



III. RESULTS 

Due to the strong in-chain coupling, any canting of 
neighboring spins away from a perfect antifcrromagnetic 
alignment within each chain will be small and unde- 
tectable in our measurements. In analyzing the data, 
we could therefore safely rely on an approximation in 
which within each chain all spins are collinear. Inter- 
chain coupling is considerably weaker than the in-chain 
one, and is comparable in strength to the Zeeman energy 
associated with experimental fields. For this reason the 
relative orientation of staggered moments in neighboring 
chains could not be considered fixed, as it may actually 
vary with field. Since there are four chains within each 
unit cell of BaCu2Si2C>7, solving the magnetic structure 



B. Spin structure at H = 4 T: Phase II 

The magnetic structure at H = 4 T was solved us- 
ing a reverse Monte Carlo analysis, which unambiguously 
pointed to a planar state with all spins in the (a, b) plane 
and a canting by an angle 4> « 80°, as shown in Fig. ^ 
on the left. This structure is a slight distortion of a co- 
linear one, with all spins parallel to the the b axis, and 
the same relative alignment of spins in adjacent chains 
as at H — 0. A least-squares refinement of the canting 
angle and the ordered moment then yielded <fi = 78(3)°, 
too = 0.18(2) /ie, with a residual x 2 = 2.4. A compar- 
ison between measured and calculated magnetic intensi- 
ties is given in Table O. Assuming <fi — 90° (6-axis colin- 
ear state), one can still get a reasonably good agreement 
with experiment, with \ 2 — 2.8. Note that the magnetic 



3 



TABLE I: Measured magnetic intensities at H = in compar- 
ison to those calculated for the c-axis colinear state visualized 



in Fig. | 
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TABLE II: Measured magnetic intensities at H — 4 T in 
comparison to those calculated for the canted state visualized 
in Fig. |(left). 
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the b axis and nearest-neighbor spins in the a direction 
aligned parallel to each other) is totally inconsistent with 
experiment, corresponding to x 2 = 158. 

IV. FIELD DEPENDENCE. 



structure of Phase II proposed in Ref. [I] (all spins parallel 
or antiparallcl to the a-axis) is totally inconsistent with 
the present experiment, corresponding to x 2 = 64, even 
when mo is optimized to best fit the data. 



C. Spin structure at H = 5 T: Phase III. 

Fewer data points were collected at H — 5 T than for 
the other two field values (see Table III). The best agree- 
ment with the data was obtained assuming a collincar 
structure with all spins pointing along the a axis (Fig. ||, 
right), and the same relative alignment of spins in ad- 
jacent chains as at H = 0. The refined value for the 
ordered moment is mo = 0.17(2) /ib, and for this model 



X — 4.5. Again we note that the spin arrangement pro- 
posed in Ref. || for the high-field phase (all spins along 



The measured field dependencies of (411), (031) and 
(051) reflections are shown in Fig. |j. The data for Phases 
II and III were analyzed using the established models. 
The value of magnetic moment mo, as well as the canting 
angle (j> for Phase II were refined to best-fit the three 
intensities at each field. The results of the fit are shown 
in solid lines in Fig. [| and the field dependence of the fit 
parameters is plotted in Fig. |]. 

To obtain a good fit to the field dependencies measured 
in Phase I, an additional parameter had to be introduced 
to account for the large change in intensities just below 
the phase transition at H c \. By trial and error, it was 
found that the best way to reproduce this effect is to as- 
sume that the spin structure, as a whole, tilts towards 
the b axis (Fig. |J). In this model, when the direction 
of staggered moment passes through the (031) direction, 
the corresponding Bragg intensity shows a characteris- 
tic dip (arrow in Fig. ||a). This behavior is a result of 
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TABLE III: Measured magnetic intensities at H — 5 T in 
comparison to those calculated for the a-axis colinear state 
visualized in Fig. ^ (right). 
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FIG. 4: Symbols: Measured field dependence of ordered 
moment (a), tilt angle a in Phase I (b) and canting angle cj> 
in Phase II (c). The solid lines are a guide for the eye. 



the neutron polarization intensity factor going to zero 
for scattering vectors that are parallel to ordered mag- 
netic moment. A smaller dip, at a slightly larger field 
is seen in the (051) peak as well (arrow in Fig. ||b). By 
varying too and the tilt angle a, a good fit to the data 
collected at < H < H c \ could be obtained. The result 
of the fit is shown in solid lines in Fig. ||, and the field de- 
pendence of the fit parameters is plotted in Fig. ^. Note 
that the angle a is shown as negative at small fields. In 
fact, the tilt a probably remains close to zero and posi- 
tive, except immediately before the transition, where it 
rapidly increases to 90° . The negative values obtained in 
the fit are an artifact of using only three measured inten- 
sities in the analysis. For small a the experiment is fairly 
insensitive to this parameter, due to a cos 2 a term in the 
intrinsic polarization factor for the scattering cross sec- 
tion and the fact that all data were collected close to the 
(001) reciprocal-space plane. However, as a approaches 
90°, the polarization sensitivity is greatly enhanced. 



FIG. 3: Measured field dependencies of several magnetic 
Bragg reflections in BaCu2Si207 (symbols). The solid lines 
are a fit to the data at each field using the models described 
in the text. The arrows indicate intensity dips just below the 
phase transition at H c \. 



V. DISCUSSION 

The data analysis described above was based on the 
known room-temperature crystal structure. The use of 
room-temperature structure factors to determine spin 
orientations at low temperatures may lead to systematic 
errors. Additional errors may be associated with measur- 
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FIG. 5: Tilting of the spin structure in BaCu2Si2 07 in mag- 
netic fields H < H c i — 2 T, applied along the c axis (Phase 
I). The spins are largely confined to the (b, c) crystallographic 
plane. 



ing magnetic intensities in the temperature-differential 
mode. Nevertheless, the results allow us to establish the 
main features of the three phases: 1) Nearest-neighbor 
spins within each chain remain almost anti-parallel to 
each other at all times. 2) The relative alignment of 
adjacent spins from neighboring chains does not deviate 
significantly from that at H = 0, with ferromagnetic spin 
alignment along a and antiferromagnetic along b. Even 
in Phase II the canting is not very dramatic and vanishes 
just before the transition at H C 2. 3) The spins are ori- 
ented roughly along the c, b and a axis in phases I, II and 
III, respectively. 4) The ordered moment is only slightly 
field-dependent and remains close to 0.15/ib- 

In weakly coupled quantum spin chains the ordered 
moment mo is proportional to y/\ J±J/J, where is 
the effective mean-field inter-chain coupling strength!] At 
H = the structure is such that Jj_| = Ji, + 2J 3 + | J a |.El 
The fact that mo remains practically constant when go- 
ing through the phase transitions, means that in Phases 
II and III, J± does not change much and nearest-neighbor 
spins remain close to being parallel to one another along 
the a axis and antiparallel along the b and c directions. 
This can be seen as an important self-consistency check 
for our determination of the spin orientations in the three 
phases. The overall gradual increase of the ordered mo- 
ment as a function of H seen in Fig. ^a is to be expected, 
and is due to the suppression of ID quantum spin fluc- 
tuations by the symmetry-breaking external field. The 
slight suppression of mo observed in Phase II (Fig. ||a) is 
a result of the canted geometry, in which Jj_ is reduced 
compared to that in Phase I. In the presence of cant- 
ing (0 < 90°) the effective inter-chain coupling is given 
by \J±\ — 2J3 — (Jft + J a \) cos(2(/>). In particular, just 
above H c \, 4> = 60°, and the drop in mo expected on go- 
ing through the first phase transition due to the abrupt 
change in |Jj_| is about 23%. This value is in reasonable 
agreement with the observed « 15% effect. 

Given the available data, the spin re-orientation at H c \ 
can be described as a more or less conventional spin-flop 
transition, that one expects for a classical antiferromag- 
net in a magnetic field applied along the easy axis. In 
zero field, the easy-c-axis anisotropy energy is minimized 



when all spins are perpendicular to the (a, b) plane. To 
take advantage of the Zeeman energy in applied fields, 
the spins have to align themselves perpendicular to c, to 
allow a canting of the rigid antiferromagnetic structure 
in the direction of H. Dzyaloshinskii-Moriya interactions 
and inter-chain coupling may modify the transition field 
and produce the small canting found in Phase II, but do 
not affect the basic physics of the effect. The fact that 
the transition appears to have a precursor, the spin struc- 
ture tilting by an angle a at H < H c \, may have several 
explanations. One possibility is that this phenomenon 
is result of an imperfect alignment of the external mag- 
netic field relative to the c axis of the crystal. Alterna- 
tively, it can be an intrinsic effect, due to the presence 
of Dzyaloshinskii-Moriya interaction, or to the existing 
structural canting of the local anisotropy axes of Cu 2+ 
relative to the c axis. To fully resolve this problem, a 
careful determination of the magnetic structure just be- 
low and just above H c \ will be required. 

The second phase transition is much more unusual, 
since it involves a spin rotation in the plane perpendic- 
ular to H. It is clear that such behavior requires the 
presence of off-diagonal exchange interactions in the sys- 
tem. As recently observed in K2V3O8, for example, a 
spin rotation around the field direction can be caused 
by a competition between DM interactions and magnetic 
anisotropy.Ll This type .of behavior is known since ear- 
lier studies of hematitefl For BaCu2Si2 0y though, a to- 
tally different explanation was proposed in Ref. o, and 
the transition was attributed to a competition between 
DM coupling in the chains and isotropic inter-chain in- 
teractions, a mechanism similar to that previously dis- 
cussed in relpfance to the anomalous spin-flop behavior 
of La 2 Cu0 4 B 

The key argument of Ref. [| is that the main com- 
ponents of the Dzyaloshinskii vector for the nearest- 
neighbor bond within the chains lie in the (a, b) plane, 
and alternate sign from one bond to the next. In 
this case, DM interactions produce a weak-ferromagnetic 
canting of the spins within each chain even in zero applied 
field. In the presence of an external field, additional cant- 
ing is due to Zeeman energy. The free energy of the sys- 
tem is minimized when the two canting effects are in the 
same direction. For this to be the case, nearest-neighbor 
spins along the a axis have to be almost antiparallel to 
each other, since Dzyaloshinskii vectors in the adjacent 
chains are antiparallel as well. This preferred antiparal- 
lel orientation is in competition with ferromagnetic inter- 
chain coupling along the a axis. In the model proposed 
in Ref. ^, the frustration is resolved in Phase II by hav- 
ing all spins aligned almost parallel to the Dzyaloshinskii 
vector (roughly along the a axis) , which practically elim- 
inates the Dzyaloshinskii energy altogether, and the fer- 
romagnetic inter-chain coupling energy along the a axis 
is minimized. In Phase III the model predicts that the 
combination of DM and Zeeman energies wins over inter- 
chain interactions: the spins rotate to be almost parallel 
to the b axis and cant in the direction of the applied field 
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to minimize both energies. At the same time, nearest- 
neighbors along the a axis become almost antiparallel, 
despite a ferromagnetic coupling in this direction. 

The results presented above clearly demonstrate that 
the actual spin arrangements in BaCu2Si20y in Phases II 
and III are totally different from those previously conjec- 
tured. In Ref. |^ it was emphasized that the Dzyaloshin- 
skii vectors not being strictly parallel to the a axis, and 
additional inter-chain interactions along the (0, 1, 0) and 
(1,1,0) directions being present, the critical field and the 
details of the spin structures may deviate from those pre- 
dicted by the simplified model. However, even in the 
more general case, if the H C 2 transition in BaCu2Si207 
was indeed driven by a competition between DM interac- 
tions and isotropic inter-chain coupling, the energy of the 
latter would have to change abruptly upon going through 
the transition. The present data clearly show that this is 
not the case, since at H C 2 the spin structure rotates as a 
whole, and the relative spin orientations, and with them 
the inter-chain exchange energies, remain practically un- 
changed. It appears that inter-chain coupling dominates 
over DM interactions in all three phases. Because of the 
canting in Phase II, the energy of inter-chain interac- 
tions is actually better minimized in Phase III, rather 
than Phase II. This sequence is opposite to that emerg- 
ing from the model proposed in Ref. || 

We believe that at the present stage it is unwise to 
attempt a quantitative explanation of the phase transi- 
tions in BaCu2Si2C>7 based on some specific spin Hamil- 
tonian. Too many parameters are involved to make such 
an analysis unambiguous. The energy scale of inter-chain 
interactions and c-axis anisotropy are indeed similar to 
that defined by transition fields. However, the Hamil- 
tonian is expected to have other terms on the same 
energy scale, and these can not be dismissed. In par- 
ticular, measurements of the spin wave spectrum con- 



tain evidence of an in-plane anisotropy of the magnitude 
« 0.15 meV. This anisotropy is most likely due to the 
so-called Kaplan Schekhtman Entin-Wohlman Aharony 
(KSEA) interactions El The KSEA term is a companion 
to DM interactions, and is an easy axis parallel to the 
Dzyaloshinskii vector. KSEA interactions are believed 
to be the driving, force in the spin reorientation transi- 
tion in K2V3O8EI The similarity of the behavior of the 
latter compound with that found in BaCu2Si207 may in- 
dicate that anisotropy effects, and KSEA interactions in 
particular, could be responsible for the H C 2 transition in 
the silicate as well. 



VI. CONCLUSION 

In summary, the actual mechanism of the exotic two- 
stage spin flop transition in BaCu2Si2C>7 remains a mys- 
tery. What is clear though, is that DM interactions 
and inter-chain exchange coupling alone can not suffi- 
cient to explain phenomenon. It appears more likely that 
both transitions in BaCu2Si207 are caused by a compe- 
tition between DM interactions and diagonal magnetic 
anisotropy effects, of which more needs to be learned. 
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